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CONVERSION  FACTORS,  NON-SI  TO  SI  (METRIC) 
UNITS  OF  MEASUREMENT 


Non-SI  units  of  measurement  used  in  this  report  can  be  converted  to  SI 
(metric)  units  as  follows: 


Multiply _ 

degrees  (angle) 
cycles  per  second 
cycles  per  second 

feet 

inches 

acceleration  of 
gravity  (standard) 


gal 

feet/second/second 

pounds 

tons 


_ By _ 

0.01745329 

1.0 

6.28318531 

0.3048 

2.54 

980.665 

32.174 

386.086 

1.0 

30.4838 

4.4822 

8.896 


To  Obtain 
radians 
hertz 

radians  per 
second 

meters 

centimeters 

centimeters/ 

second/second 

feet/second/ 

second 

inches/second/ 

second 

centimeters/ 

second/second 

centimeters/ 

second/second 

newtons 

kilonewtons 
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INTRODUCTION  TO  THE  COMPUTATION  OF  RESPONSE  SPECTRUM 
FOR  EARTHQUAKE  LOADING 


PART  I :  INTRODUCTION 

1.  This  paper  presents  an  introduction  to  the  computation  of  a  response 
spectrum  for  earthquake  loading.  A  response  spectrum  is  a  graphical  relation¬ 
ship  of  maximum  values  of  acceleration,  velocity,  and/or  displacement  response 
of  an  infinite  series  of  elastic  single  degree  of  freedom  (SDOF)  systems 
subjected  to  a  time  dependent  excitation.  To  accomplish  this  task,  the  formu¬ 
lation  and  solution  of  the  equation  of  motion  for  a  damped  SDOF  system 
subjected  to  dynamic  excitation  is  reviewed  prior  to  discussion  of  the 
response  spectrum. 

2.  Examples  of  idealized  SDOF  systems  are  shown  in  Figure  1.  The 
dynamic  excitation  may  be  due  to  the  forcing  function  P(t)  (Figure  lc)  acting 
on  the  mass  or  to  ground  shaking,  typically  expressed  in  terms  of  a  ground 
acceleration  time  history,  as  shown  in  Figure  2.  The  dynamic  response  of 
damped  SDOF  systems  is  described  by  the  variations  of  displacement,  velocity, 
and  acceleration  of  the  mass  with  time.  A  plot  of  the  maximum  values  of 
acceleration,  velocity,  and/or  displacement  of  an  infinite  series  of  SDOF  sys¬ 
tems  versus  undamped  natural  period  is  called  a  response  spectrum. 

3.  The  response  spectrum  is  the  cornerstone  of  modern  earthquake  engi¬ 
neering  and  structural  dynamics.  It  is  used  to  calculate  the  dynamic  response 
of  multi-degree  of  freedom  (MDOF)  semidiscrete  structural  models  of  buildings 
and  hydraulic  structures  (i.e.  dams,  locks,  and  intake  towers)  as  well  as  for 
the  evaluation  of  frequency  content  of  recorded  accelerograms.  The  frequency 
content  of  accelerograms  is  of  importance  for  selecting  the  ground  motion(s) 
to  be  used  as  the  forcing  functions  in  the  dynamic  analysis. 

4.  The  damped  SDOF  system  and  the  dynamic  equilibrium  equation  for  the 
system  are  introduced  in  Part  II.  The  derivation  of  the  equation  of  motion 
for  a  damped  SDOF  system  subjected  to  a  time  dependent  force  history  is 
described  in  Part  III.  Part  IV  derives  the  equation  of  motion  for  the  case  of 
a  damped  SDOF  system  subjected  to  earthquake  shaking.  This  part  also 
describes  how  the  earthquake  shaking  problem  may  be  reformulated  as  an  equiva¬ 
lent  force  history  problem. 
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(a)  Idealized  rigid  mass  SDOF  systems 
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(b)  Idealized  damped  SDOF  systems 


P 
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(c)  Variation  in  force  with  time 


Figure  1.  Idealized  single  degree  of  freedom  (SDOF)  systems 


5.  The  fundamentals  of  the  solution  of  the  equation  of  motion  using 
numerical  procedures  are  introduced  in  Part  V.  Because  of  the  irregular 
nature  of  the  acceleration  time  histories  that  have  been  recorded  during 
earthquakes,  numerical  methods  are  used  to  compute  the  response  of  SDOF  sys¬ 
tems  for  developing  response  spectra. 

6.  Part  VI  describes  the  construction  of  response  spectra  using  the 
accelerations,  velocities,  and  displacements  obtained  from  the  numerical 
integration  of  the  equations  of  motion  for  an  infinite  series  of  SDOF  systems 
subjected  to  earthquake  shaking.  The  terms  associated  with  response  spectra 
are  also  defined  in  Part  VI. 
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PART  II:  THE  DYNAMIC  EQUILIBRIUM  EQUATION  FOR  A  DAMPED  SDOF  SYSTEM 


7.  A  damped  SDOF  system  consists  of  a  rigid  mass  m,  a  linear  spring  of 
stiffness  k,  and  a  viscous  damper  with  damping  coefficient  c,  as  shown  in  Fig¬ 
ure  2.  The  viscous  damper  represents  the  energy  absorbing  component  of  the 
system.  Each  rigid  mass  is  constrained  to  move  in  the  horizontal  x  direction. 
The  SDOF  system  in  Figure  2a  is  subjected  to  a  time  dependent  horizontal  force 
history  P(t),  while  the  SDOF  system  in  Figure  2b  is  shaken  by  a  time  dependent 
horizontal  ground  acceleration,  xground(t).  The  double  overdot  represents  the 
second  derivative  of  displacement  x  with  respect  to  time.  The  dynamic 
response  of  either  system  at  any  time  t  is  governed  by  the  relationship 

E  F(  t)  =  m  xtotal  ( t)  (1) 


where 

xtotai(c)  “  the  acceleration  of  the  mass  m. 

8.  In  earthquake  engineering  problems,  the  displacements,  velocities 
and  accelerations  of  the  rigid  mass  differ  from  those  of  the  ground  (Fig¬ 
ure  2b) .  It  is  convenient  to  describe  the  motion  of  the  rigid  mass  in  terms 
of  the  relative  displacement  of  the  mass  with  respect  to  the  ground  as 

*(t>  =xcocal(t)  -  xgTound(t)  (2) 

where 

x(t)  -  the  relative  displacement 

xtotai ( t )  ”  the  displacement  of  the  mass  from  its  at-rest  position 
xground(t)  -  the  displacement  of  the  ground 
The  relative  velocity  and  relative  acceleration  of  the  mass  m  at  time  t  are 
obtained  by  differentiation  of  Equation  2 

*<t)  =  xcocal{t)  -  xgiound(t)  (3) 

and 

X{t)  =  Xtotal{t)  -  Xground{t)  W 

The  solutions  to  the  two  dynamic  problems  shown  in  Figure  2  are  described  in 
Part  III  for  an  applied  dynamic  force  and  in  part  IV  for  an  earthquake 
excitation . 
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PART  III:  ORDINARY  DIFFERENTIAL  EQUATION  FOR  A  DAMPED  SDOF  SYSTEM  - 

FORCED  VIBRATIONS 


9.  This  part  describes  the  derivation  of  the  ordinary  differential 
equation  of  second  order  that  governs  the  dynamic  response  of  the  damped  SDOF 
system  shown  in  Figure  3  to  a  time  dependent  external  force  P(t).  This  simple 
mechanical  system  serves  as  an  introduction  to  the  fundamentals  of  the 
dynamics  of  a  SDOF  system  and  as  background  information  for  the  dynamic  prob¬ 
lem  involving  a  SDOF  system  undergoing  a  time  dependent  ground  acceleration  as 
described  in  Part  IV. 

10.  Assume  that  at  time  t,  P(t)  acts  in  the  positive  x  direction  (to 
the  right)  and  that  the  acceleration,  velocity,  and  displacement  of  the  mass 
are  positive.  For  this  problem,  the  ground  is  at-rest  (i.e.  xground  -  0), 
therefore  xtotal(t)  is  equal  to  x(t).  The  movement  of  the  mass  from  its  at- 
rest  position  to  the  right  stretches  the  spring,  resulting  in  a  restoring 
force  fk(t).  For  a  linear  spring 

fk(t)  =  k  xit)  (5) 


where 

k  -  the  spring  stiffness 
x(t)  -  the  displacement  of  the  mass  at  time  t 
Similarly,  the  restoring  force  fc(t)  contributed  by  the  viscous  damper  at  time 
t  is  given  by 

fcit)  =  c  x(t)  (6) 

where 

c  -  the  damping  coefficient 
x(t)  -  the  velocity  of  the  mass  at  time  t 
Both  fk(t)  and  fc(t)  act  in  the  negative  x  direction  (to  the  left)  at  time  t, 
tending  to  restore  the  mass  to  its  at-rest  position  (refer  to  Figure  3).  The 
sum  of  forces  acting  on  the  mass,  the  left-hand-side  of  Equation  1,  is 

L  Fit )  =  Pit)  -  fkit)  -  fcit)  (7) 

Substituting  Equations  5,  6,  and  7  into  Equation  1  with  x(t)  -  xtotal(t)  and 
rearranging  results  in 
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pm 


FIXED  BASE 


2  Forces  ■  m”(t) 


DIRECTION  OF 
ACCELERATION 
OF  MASS  m 
AT  TIME  t 


mx(t)  ♦  cx(t)  ♦  kx(t)  •  P(t) 


Figure  3.  Forces  acting  on  a  linear  SDOF  system  at  time  t,  external 

force  P(t)  applied 


rnx(C)  *  c  x(t)  *  k  x(t)  =  Pit)  (8) 

The  first  term  in  Equation  8  represents  the  inertial  force  f^t),  associated 
with  the  mass  m  undergoing  an  acceleration  x(t),  which  acts  opposite  to  the 
direction  of  the  acceleration  of  the  mass  m  (refer  to  Figure  4) . 


Free  Vibration 


11.  Free  vibration  results  from  the  application  of  an  initial  displace¬ 
ment  and/or  velocity  with  no  external  forcing  function  acting  on  the  system 
(P(t)  -  0). 


12 


Ptt) 


DIRECTION  OF 
ACCELERATION 
OF  MASS  m 
AT  TIME  t 


FIXED  BASE 


fk(U 

fc(t) 


f i (t)  ♦  fc(t)  *  fk(t)  ■  P(t) 
where  the  inertiol  force  f;  is  given  by- 
f;  (t)  -  mx(t) 

Figure  4.  Inertial  force  acting  opposite  to  the 
acceleration  of  mass  m  at  time  t,  external  force 

P(t)  applied 

Undamped  Free  Vibration 


12.  The  equation  of  motion  (Equation  8)  for  an  undamped  system  (c  -  0) 
in  free  vibration  is 


m  x(t)  *  k  x(  t)  =  0 


(9) 


and  its  solution  is 


x(t) 


x„  cos  ut  ♦  —  sm  tot 
0  to 


(10) 


where 

xQ  -  the  initial  displacement 
xQ  -  the  initial  velocity 
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and 


(a) 


(11) 


where 

u  -  constant  circular  frequency  (rad/sec). 

Equation  10  is  obtained  using  standard  solution  procedures  for  linear  differ¬ 
ential  equations,  such  as  the  method  of  undetermined  coefficients  (Sec¬ 
tion  2.12  in  Kreyszig  1972).  The  motion  described  by  Equation  10  and  shown  in 
Figure  5a  is  cyclic  with  a  constant  maximum  amplitude 


l*(t)  I  max 


(12) 


and  constant  circular  frequency  u.  The  cyclic  nature  of  vibration  may  also  be 
expressed  in  terms  of  the  natural  period  of  vibration  ,  T  (sec)  ,  of  the 
undamped  SDOF  system  as 

T  =  m  (id 

0i 


or  the  natural  cyclic  frequency  of  vibration  f  (cycles/sec  or  Hz)  where 


f  = 


oj 

2  n 


(14) 


13.  Figure  6  shows  the  computed  response  of  a  SDOF  system  with  k  - 
1  N/m  and  mass  m  -  1  kg,  subjected  to  an  initial  displacement  x„  -  1  cm  (x„  - 
0).*  The  circular  frequency  is  equal  to  1  radian/sec  (by  equation  11)  and  the 
undamped  natural  period  of  vibration  is  equal  to  2  n  seconds.  Plots  of  the 
displacement,  velocity  and  acceleration  of  the  mass  for  the  first  cycle  of 
harmonic  response  are  shown  in  Figure  6,  as  well  as  the  variation  with  time  of 
the  spring  force  and  inertial  force  acting  on  the  rigid  mass.  Note  that  the 
inertial  force  acts  opposite  to  the  direction  of  the  acceleration  vector  for 
the  rigid  mass. 


*  One  newton,  N,  is  the  force  which  gives  an  acceleration  of  1  m/sec2  to  a 
mass  of  1  kg. 
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DISPLACEMENT 


PERIOO  OF  VIBRATION:  T 


2* 


♦x 

0 

-x 


time 

(sec) 


(a)  Undamped  free  vibration 


x  . 


x  ♦  20m  i  ♦  oj2  x  -  0 


PERIOO  OF  VIBRATION:  TD  *  t/ThS7 


(b)  Damped  free  vibration  (damped  at  a  ratio  to  critical 

damping  equal  to  0) 

Figure  5.  Free  vibration  response  of  damped  and  undamped  SDOF 

systems 
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f  k  =  kx  J*  2  0 

-0.01 

0.01 

f;  =  mx*  ^.-z  0 

-0.01 


Figure  6.  Undamped  free  vibration  responses  of  a  SDOF  system 
given  an  initial  displacement  of  1  cm 
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Damped  Free  Vibration 


14.  The  equation  of  motion  for  damped  free  vibration  may  be  written  as 


St(t)  +  2  p  w  x(t)  *  x(t)  =  0 


(15) 


where 


P  = 


c 

2  mu 


(16) 


Equation  15  was  derived  from  Equation  8  by  (1)  dividing  each  term  by  the  mass 
m,  (2)  introducing  the  constant  /}  for  the  damping  force  term,  and  (3)  P(t)  set 
equal  to  zero.  The  solution  of  Equation  15  is,  for  0  <  1, 


x(t) 


COS  <1>DC  * 


X0  ♦  p<dx0 
<*>0 


sin  v>Dt 


\ 


(17) 


The  motion  described  by  Equation  17  for  0  <  1,  shown  in  Figure  5b,  is  periodic 
with  an  exponentially  decaying  amplitude  and  damped  circular  frequency 

o)D  =  u  -  P*  (18) 


or  period 


Td 


T 


(19) 


where  w  and  T  are  the  circular  frequency  and  period  of  the  undamped  system, 
respectively. 

15.  When  p  -  l  the  period  of  the  damped  system  becomes  infinite  indi¬ 
cating  that  the  motion  is  no  longer  cyclic  but  after  an  initial  maximum  dis¬ 
placement,  decays  exponentially  without  subsequent  reversal  of  direction.  A 
system  with  f)  -  1  or  with  damping 

*  =  cCIltic.j  =  2  m  o»  (20 


is  said  to  be  "critically"  damped.  Hence 
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p » 


c 

critical 


C 

2  m  <ji 


(21) 


is  called  the  "damping  ratio"  or  the  "fraction  of  critical  damping." 

Typically  for  structures  0  is  less  than  0.1  and  for  material  damping  in  soils 
is  less  than  0.25.  may  also  be  expressed  as  percentage.  For  typical  struc¬ 
tural  systems  the  periods  (or  frequency)  of  the  damped  and  undamped  systems 
are  approximately  equal,  i.e. 

T  <=  Td  (22) 


Forced  Vibration  with  Dynamic  Force  Applied  to  Mass 

! 

16.  The  equation  of  motion  (Equation  8)  for  the  SDOF  system  (Figure  3) 
may  be  written  as 

*(t)  ♦  2  u  *(t)  +  x(t)  =  (23)  | 

m  i 


In  general  there  are  three  approaches  to  the  solution  of  Equation  23:  closed 
form  solutions,  Duhamel’s  Integral,  and  numerical  methods. 

Closed  Form  Solutions 

17.  For  a  simple  harmonic  force  history,  such  as  P(t)  -  Constant- sin 
[“driv*'  t]  in  Equations  8  and  23,  closed  form  solutions  are  available  in  numer¬ 
ous  textbooks  on  both  mechanical  vibrations  and  structural  dynamics.  This 
procedure  is  not  practical  for  earthquake  engineering  problems  involving 
irregular  force  time  histories. 

Duhamel's  Integral 

18 .  A  second  procedure  used  to  solve  for  the  dynamic  response  of  a  SDOF 
system  involves  the  representation  of  the  load  time  history  P(t)  as  a  series 


18 


of  impulse  loadings  P(r)  that  are  applied  to  the  SDOF  system  for  infinitesimal 
time  intervals  dr,  as  shown  in  Figure  7. 

x  . 


h*i 

k 

1 - 

\ 

1 

HASS,  m 

i 

\ 

1 

\pm 

i 

i 

dr  ■  TIME  INTERVAL 

Figure  7.  Load  time  history  P(r)  as  impulse  loading  to  undamped  SDOF 

system 

19.  Figures  8a  through  8c  outline  the  derivation  of  the  incremental 
displacement  dx  of  an  undamped  SDOF  system  at  time  t  to  a  single  force  pulse 
P(r)-dr.  Figure  8d  gives  the  resulting  integral  solution,  known  as  the 
Duhamel's  integral,  for  the  displacement  of  the  undamped  SDOF  system  at  time  t 
to  the  s recession  of  force  pulses  applied  between  time  -  0  and  time  -  t. 

20.  Duhamel's  integral  for  a  damped  SDOF  system  is 

x  { t)  -  — l—  fC  P(  x)  sin  [®„(t  -  r)l  dx  (24) 

m  J 


where 

wt,  -  the  damped  circular  frequency  of  vibration  (Equation  18)  , 
u  -  the  undamped  circular  frequency  of  vibration  (Equation  11) 

P  -  the  fraction  of  critical  damping  (Equation  21). 

21.  For  a  select  few  load  time  histories  P(t),  the  integral  may  be 
evaluated  directly.  Irregular  force  time  histories  require  numerical  solu¬ 
tions  to  be  used  to  evaluate  Duhamel's  integral.  These  procedures  are 
described  in  numerous  textbooks  on  structural  dynamics. 

Numerical  Methods 

22.  For  most  earthquake  engineering  problems,  numerical  methods  are 
used  to  solve  a  form  of  Equation  8  or  Equation  23  due  to  the  irregular  form  of 
force  time  histories.  Numerical  methods  will  be  discussed  in  Part  V. 
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Jr  , 


1 

1 

1 

\ 

!  pm 

HASS.m 

\ 

t 

i 

i 

l)'  ]  ir 

ii  i 

FIXED  BASE 

IMPULSE  -  CHANGE  IN  MOMENTUM 

P(Ddr  -  m(dx)  or  dx  -  — —  P(r)dr 

m 

dx  -  INCREMENTAL  VELOCITY  -  x0(D 


P 


P(TH 


IMPULSE  *  PiTi'dT 


d  r 


dr  ) 


a  -t) 

TIME  INTERVAL 


(a)  Incremental  velocity  dx  at  time  r 


mB  ♦  kx  -  0 

INITIAL  CONDITIONS:  x(t-O)  •  0 
x(t-O)  -  x0 


(b)  Exact  solution  for  free  vibration  response 

(undamped) 


Figure  8.  Duhamel's  integral  for  an  undamped  SDOF 

system 


-  time 
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Introducing  dx  from  <o)  into  x0  of  (b),  the  incrementol  displocement  is  given  by 
dxtt)  -  -1-  P(T)<jrj  sin  [um  -  rTj 

(c)  Incremental  displacement  at  time  t  due  to  a  single  pulse 

at  time  r 


xtt)  •  dxtt)  dT 


X<°  ’  wjf  PlT)  s!n[^(t  *r^j 

(d)  Duhamel's  integral 
Figure  8.  (Continued). 
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PART  IV:  ORDINARY  DIFFERENTIAL  EQUATION  FOR  A  DAMPED  SDOF  SYSTEM  - 

GROUND  ACCELERATION 

23.  This  part  describes  the  derivation  of  the  ordinary  differential 
equation  of  second  order  that  governs  the  dynamic  response  of  the  damped  SDOF 
system  shown  in  Figure  2b  shaken  by  a  horizontal  ground  acceleration 
xground(t)-  Aground  varies  in  both  magnitude  and  direction  with  time.  Assume 
that  at  time  t,  xground(t)  is  positive  (to  the  right)  and  the  acceleration, 
velocity,  and  displacement  of  the  mass  are  positive.  The  restoring  forces  of 
the  spring  and  dashpot  shown  in  Figure  9  are  given  by  Equations  5  and  6, 
respectively.  Combination  of  Equation  1 

E  F(t)  =  mxcoei,j(t)  (1) 

Equation  5 

fk(t)  -  k  x  (t)  (5) 

and  Equation  6 

fc(t)  =  c  x(t)  (6) 

result  in 

m  keocal(t)  +cx(t)  +  k  x[t)  =0  (25) 

where  x  and  x  are  relative  velocity  and  displacement  of  the  mass  (Equations  2 
and  3)  valid  at  any  time  t.  The  first  term  in  Equation  19  represents  the 
inertial  force  f L ,  associated  with  the  mass  m  undergoing  a  total  acceleration 
xtotal(t).  This  force  vector  acts  opposite  to  the  direction  of  the  total 
acceleration  vector  of  the  mass  m,  as  shown  in  Figure  9.  Substituting  Equa¬ 
tion  4  into  Equation  25  for  xtotal  results  in  the  expression 

m  Jt{t)  *  c  x(t)  *  k  x  =  -  m  xgTOUnd{t)  (26) 

or 

*(t)  *2  Pox(t)  ♦  x(t)  =  -  xground{t)  (27) 


AT-REST 


k 

c_^ 

i 

UASS,  m 

.  C 

L? 


ground 


time 


fk(t) 

fe(t) 


f;(t)  ♦  fc(t)  -  fk(t>  -  0 
WHERE  •• 

f;(t)  -  mxt0t0|(t) 

'  mC*ground(t)  *  x(t)] 

Figure  9.  Forces  acting  on  an  SDOF  system  at  time  t, 
ground  acceleration  xground 

where 

w  -  the  undamped  circular  frequency  (Equation  11) 

-  the  damping  ratio  (Equation  21) 

The  terms  associated  with  the  relative  movement  of  the  mass  are  collected  on 
the  left-hand-side  of  the  equal  sign.  Either  differential  equation 
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(Equation  25  or  26)  describes  the  dynamic  response  of  a  damped  SDOF  system 
shaken  by  a  horizontal  ground  acceleration  xground. 


Equivalent  SDOF  Problems 


24.  Comparison  of  Equation  26  with  Equation  8  shows  that  the  relation¬ 
ships  for  the  two  SDOF  systems  shown  in  Figure  2  differ  by  the  term  on  the 
right -hand- side  of  the  equal  sign  (the  force  history).  Thus,  the  problem  of  a 
damped  SDOF  system  shaken  by  a  time  varying  ground  acceleration  is  equivalent 
to  the  problem  of  a  damped  SDOF  system  resting  on  a  fixed  base  and  subjected 
to  a  time  dependent  force  P(t)  of  magnitude  -  rn-xground,  as  shown  in  Figure  10. 


STEP  1.  SOLVE 

mx(t)  ♦  cx(t)  ♦  kx(t)  ■  *mxground(t) 
STEP  2,  SOLVE 

"  x(t)  *  X  ground^ 

Figure  10.  Equivalent  dynamic  SDOF  system  problems 


Solution  of  the  Equation  of  Motion 


25.  The  total  dynamic  response  of  the  SDOF  system  problem  is  computed 
in  two  steps.  Step  1  solves  for  the  relative  response  of  the  damped  SDOF 
system  as  governed  by  the  ordinary  differential  Equation  27,  and  in  Step  2, 
the  total  response  is  equal  to  the  sum  of  the  relative  response  (from  step  1) 
plus  the  motion  of  the  ground. 
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Closed  Form  Solutions 


26.  For  simple  harmonic  ground  accelerations  (e.g.  xground  -  Con¬ 
stant-  sin  [wdrive-  t]  )  closed  form  solutions  to  Equation  27  are  available  in 
numerous  textbooks  on  both  mechanical  vibrations  and  structural  dynamics. 

This  procedure  is  impractical  for  earthquake  engineering  problems  due  to  the 
irregular  nature  of  ground  acceleration  time  histories. 

Duhamel's  Integral 

27.  A  second  procedure  used  to  solve  for  the  relative  displacement  of 
the  SDOF  system  involves  the  representation  of  the  load  time  history  P(t)  -  - 
mXground ( t )  as  a  series  of  impulse  loadings  P(t)  applied  to  the  SDOF  system 
for  infinitesimal  time  intervals  dr  (Figure  7) .  By  introducing  P(t)  -  - 
mXground ( t )  into  Equation  24,  Duhamel's  integral  for  a  damped  SDOF  system  is 

X(t)  =  -  f  xground(x)  sin  [o>D(t  -  t)  ]  dx  (28) 

u)  n  * 

D  0 


where 

-  the  damped  angular  frequency  of  vibration  (Equation  18) 

«  -  the  undamped  angular  frequency  of  vibration  (Equation  11) 

P  -  the  fraction  of  critical  damping  (Equation  16) 

The  irregular  forms  of  acceleration  time  histories  require  numerical  solutions 
to  be  used  to  evaluate  Duhamel's  integral. 

Numerical  Methods 


28.  In  usual  applications  to  earthquake  engineering  problems,  numerical 
methods  are  used  to  solve  Equation  27  or  Equation  28  for  the  relative  dis¬ 
placement  of  the  SDOF  mass  because  of  the  irregular  nature  of  ground  accelera¬ 
tion  time  histories  (discussed  in  Part  V) . 
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PART  V:  SOLUTION  OF  DYNAMIC  EQUILIBRIUM  EQUATIONS  FOR  A  DAMPED 
SDOF  SYSTEM  USING  NUMERICAL  METHODS 


29.  This  part  introduces  the  fundamentals  of  numerical  methods  used  to 
solve  for  the  accelerations,  velocities,  and  displacements  of  a  damped  SDOF 
system  due  to  a  time  dependent  loading.  In  general,  there  are  two  categories 
of  numerical  methods  used  for  solving  the  dynamic  equilibrium  equation: 
direct  integration  methods  and  frequency -domain  methods.  This  part  discusses 
direct  integration  methods  only.  The  reader  is  referred  to  books  on  struc¬ 
tural  dynamics  for  a  description  of  frequency -domain  methods. 


Direct  Integration  Methods 


30.  Direct  integration  methods  are  used  to  solve  for  the  response  of 
the  SDOF  (and  MDOF  semidiscrete  structural  models)  by  direct  integration  of 
the  dynamic  equilibrium  equations  at  closely  spaced,  discrete  time  intervals 
throughout  the  time  of  shaking  using  a  numerical  step-by-step  procedure  of 
analysis.  The  term  "direct"  means  that  prior  to  numerical  integration,  there 
is  no  transformation  of  the  equations  into  a  different  form,  as  is  done  in  a 
frequency-domain  analysis.  Table  1  lists  some  of  the  step-by-step  algorithms 
used  in  structural  dynamics  for  both  SDOF  systems  and  MDOF  semidiscrete  struc¬ 
tural  models  and  in  the  characterization  of  ground  motions  for  earthquake 
engineering  problems. 

31.  Direct  integration  methods  are  based  on  two  concepts.  First,  the 
equation  of  motion  (Equation  27)  is  satisfied  at  discrete  points  in  time  (i.e. 
t,  t  +  At,  ...)  during  earthquake  shaking,  and  second,  the  forms  of  the  varia¬ 
tion  in  displacements,  velocities,  and  accelerations  within  each  time  inter¬ 
val,  At,  are  assumed.  Direct  integration  time  methods  are  classified  as 
either  explicit  integration  methods  or  implicit  integration  methods.  The 
explicit  integration  method  solves  for  the  unknown  values  of  xt  +  At,  xt  +  At, 
and  xt  +  At  at  each  new  time  (t  +  At)  using  the  equation  of  motion  at  time  t, 
with  the  known  values  for  xt,  xt,  and  xt  at  time  t  as  the  initial  conditions. 
The  implicit  integration  method  solves  for  the  unknown  values  of  xt  +  At,  xt  +  At 
and  xt  +  At  at  each  new  time  (t  +  At)  using  the  equation  of  motion  at  time  (t  + 
At) .  For  MDOF  systems  implicit  schemes  require  the  solution  of  a  set  of 
simultaneous  linear  equations,  whereas  explicit  schemes  involve  the  solution 


Table  1 

Step-bv-Step  Algorithms  Used  In  Structural  Dynamics 


Family  of 
Structural 
Dynamics 
Aleorithms 

Example 

of 

Algorithm 

Type 

Stability 

Condition 

Order  of 
Accuracv 

collocation 

methods 

Wilson-0 

Implicit 

unconditional 
for  0>1.366 

0(At)2 

Newmark-/} 

methods 

Average  acceleration 
(trapezoidal  rule) 

Implicit 

unconditional* 

O(At)2 

Linear  acceleration 

Implicit 

conditional 

O(At)2 

Fox-Goodwin  formula 

Implicit 

conditional 

o 

/•“N 

> 

ft 

fo 

Central  difference 

Explicit 

conditional 

O(At)2* 

Houbolt' s 
method 

Implicit 

unconditional 

O(At)2 

a -method 

Hilber -Hughes -Taylor 

Implicit 

unconditional* 

o 

> 

ft 

to 

Wood-Bossak- 

Zienkiewicz 

Implicit 

unconditional* 

O(At)2 

-method 

Hoff-Pahl 

Implicit 

unconditional* 

O(At)2 

beta-m 

method 

Katona-Zienkiewicz 

Implicit 

unconditional* 

O(At)2 

*  For  select  values  of  constants  used  in  algorithm. 


of  a  set  of  linear  equations,  each  of  which  involves  a  single  unknown.  Both 
implicit  and  explicit  step-by-step  algorithms  are  listed  in  Table  1.  Implicit 
algorithms  are  the  more  popular  of  the  two  types  of  numerical  methods  in 
earthquake  engineering  problems  because  of  the  larger  size  time  step  that  may 
be  used  in  the  analysis.  However,  implicit  methods  involve  considerable  com¬ 
putational  effort  at  each  time  step  compared  to  explicit  methods  since  the 
coefficient  matrices  for  MDOF  systems  must  be  formulated,  stored,  and  manipu¬ 
lated  using  matrix  solution  procedures. 
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Linear  Acceleration  Method 


32.  The  principles  common  to  all  step-by-step  time  integration  solu¬ 
tions  of  the  equations  of  motions  are  illustrated  for  the  Figure  10  SDOF  sys¬ 
tem  by  applying  the  linear  acceleration  algorithm  to  this  problem  for  a  single 
time  interval  At.  The  dynamic  response  of  the  mass  to  an  earthquake  time 
history  at  each  time  step  is  expressed  in  terms  of  the  values  for  the  dis¬ 
placement,  the  velocity,  and  the  acceleration  of  the  mass. 

33.  The  linear  acceleration  algorithm,  one  of  the  simpler  forms  of  the 
Newmark-/J  family  of  implicit  algorithms,  assumes  a  linear  variation  in 
acceleration  of  the  mass  from  time  t  to  time  (t  +  At)  as  shown  in  Figure  11. 
The  values  for  the  three  variables  are  known  at  time  t,  and  the  values  are 
unknown  at  time  (t  +  At). 

34.  The  assumed  linear  variation  in  acceleration  provides  one  of  the 
four  equations  used  in  the  linear  acceleration  algorithm.  The  slope  is 
expressed  in  terms  of  the  values  of  acceleration  at  time  t  and  time  (t  +  At) , 
as  listed  in  Figure  11.  Two  additional  equations  are  provided  by  twice  inte¬ 
grating  the  linear  acceleration  relationship  from  time  t  to  time  (t  +  At). 

This  results  in  a  quadratic  variation  in  velocity  of  the  mass  over  time  step 
At  and  a  cubic  variation  in  displacements  of  the  mass  over  time  At  (Fig¬ 
ure  11).  The  fourth  equation  is  given  by  the  equation  of  motion  at  time  (t  + 
At)  . 

35.  Figure  12  summarizes  the  three  steps  when  solving  for  the  dynamic 
response  of  the  mass  at  each  new  time  (t  +  At) .  The  three  relationships  shown 
in  this  figure  were  obtained  by  rearranging  the  four  relationships  listed  in 
Figure  10.  The  first  stage  of  the  analysis  involves  the  application  of  the 
step-by-step  procedure  of  analysis  during  the  time  of  earthquake  shaking, 
computing  the  time  histories  of  response  for  the  mass.  Due  to  the  nature  of 
the  formulation,  the  computed  acceleration  time  history  is  the  relative 
acceleration  of  the  mass.  The  total  acceleration  of  mass  is  equal  to  the  sum 
of  the  relative  acceleration  values  computed  at  each  time  step,  and  the  corre¬ 
sponding  ground  acceleration  values  (step  2  in  Figure  10). 

Time  Step 

36.  The  selection  of  the  size  of  the  time  step  At  to  be  used  in  the 
step-by-step  calculation  of  the  dynamic  response  of  the  SDOF  (and  of  MDOF 


Known 

Unknown 

Value 

Value 

Variable 

(time  -  t) 

(time  -  t  * 
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Figure  11.  Step-by-step  time  integration  solution  of  the  equation  of 
motion  using  the  linear  acceleration  method 
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LHS  (unknown) 


RHS  (known) 


step  l,  solve  constant ^  xc+At  =  constant ^ 
step  2,  solve  for  xt+At  =  -^.]x£+At  -  termx 

step  3,  solve  for  *t+Ac  =  -----  x£+At  -  term. 2 

,  (At)2. 

where  the  constants  are  given  by 

constant tJt<i  =  k  +  f-^-1  c  +  — - —  m 

[At]  ( A  t) 2  J 

constant ^  =  P£tA£  +  cterrr^  +  mterm2 

and 

~  tn  ■^yroun<j(t  +  A  t) 
ter/77i  -  *t  +  2  xt  +  A£  *£ 


Figure  12.  Solution  for  the  dynamic  response  of  the  mass  at  time 
(t  +  At)  using  the  linear  acceleration  algorithm 
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semidiscrete  structural  models)  is  restricted  by  stability  and  accuracy  con¬ 
siderations.  The  primary  requirement  of  a  numerical  algorithm  is  that  the 
computed  response  converge  to  the  exact  response  as  At  -*  0  (Hughes  1987).  The 
stability  and  accuracy  criteria  are  expressed  in  terms  of  a  maximum  allowable 
size  for  the  time  step,  Atcritlcal,  which  differs  among  the  various  numerical 
algorithms . 


Stability 


37.  The  stability  condition  requirements  for  numerical  algorithms  are 
categorized  as  either  unconditional  or  conditional  (Table  1) .  Bathe  and 
Wilson  (1976),  Bathe  (1982),  and  Hughes  (1987)  describe  an  integration  method 
as  unconditionally  stable  if  the  numerical  solution  for  any  initial  value 
problem  (e.g.  Figure  5a)  does  not  grow  without  bound  for  any  time  step  At, 
especially  if  the  time  step  is  large.  The  method  is  conditionally  stable  if 
the  previous  statement  is  true  only  for  those  cases  in  which  At  is  less  than 
some  critical  time  step  Atcritical.  The  stability  criterion  for  a  numerical 
algorithm  is  established  by  the  values  assigned  to  the  constants  that  are  used 
in  the  algorithm  and  the  terms  associated  with  the  structural  model  (Hughes 
and  Belytshko  1983,  Hughes  1987,  and  Dokainish  and  Subbaraj  1989).  For 
example,  numerical  stability  considerations  for  the  linear  acceleration  method 
require  that  the  time  step  At  of  a  SDOF  system  be  restricted  to  values  given 
by  the  relationship 


At  s  A tcrieical 


(29) 


where 


A  t 


critical 


=  1  =  T  & 

n  f  n  n 


(30) 


The  values  of  Aterltlcal  are  summarized  in  Hughes  and  Belytshko  (1983),  Hughes 
(1987),  and  Dokainish  and  Subbaraj  (1989)  for  other  algorithms. 

38.  The  attributes  of  a  stable  numerical  analysis  are  illustrated  using 
the  free  vibration  problem  of  an  undamped  SDOF  system  shown  in  Figure  13, 
given  an  initial  displacement  x0  and  an  initial  velocity  Xq.  With  no  damping, 
the  exact  solution  for  this  initial  value  problem  is  a  harmonic  function  that 
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DISPLACEMENT 
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- 1 
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[  INITIAL  CONDITIONS:  x(t-O)  ■  x0 

i  x(t*0)  *  5, 

EXACT  SOLUTION:  x  •  x0  cos  u>t  ♦  —  sin  cat 


where1  u> 


(a)  Stable  free  vibration  response  with  a  smaller  time 
step  than  the  critical  time  step 


(b)  Unstable  free  vibration  response  due  to  a  larger  time 
step  than  the  critical  time  step 

Figure  13.  Example  of  unstable  response  for  an  undamped 
SDOF  system  in  free  vibration 
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is  continuous  with  time  and  equal  to  the  sum  of  a  sine  wave  plus  a  cosine  wave 
of  constant  maximum  amplitude.  Figure  13a  shows  a  conditionally  stable  numer¬ 
ical  solution  that  corresponds  to  the  exact  solution.  Figure  13a  results 
contrast  with  Figure  13b  results  for  the  unstable  free  vibration  response 
computed  using  At  >  Atcritical,  in  which  the  maximum  computed  amplitude  of  the 
SDOF  system  increases  with  time. 

39.  The  value  of  Atcritlcal  for  MDOF  systems  is  governed  by  the  smallest 
period  or  the  largest  frequency  of  the  semidiscrete  MDOF  model.  Violation  of 
the  stability  criterion  by  an  even  single  mode  will  destabilize  the  numerical 
solution  of  the  dynamic  problem.  To  avoid  this  restriction  and  increase  the 
time  step  size,  unconditionally  stable  integration  methods  are  used. 

Accuracy 

40.  The  accuracy  of  a  numerical  algorithm  is  associated  with  the  rate 
of  convergence  of  the  computed  response  to  the  exact  response  as  At  -*  0 
(Hughes  1987).  For  the  dynamic  analysis  of  MDOF  semidiscrete  structural 
models  using  the  finite  element  method,  the  accuracy  of  an  algorithm  is  also 
concerned  with  the  question  of  the  computed  response  of  spurious  higher  (fre¬ 
quency)  modes  of  the  semidiscrete  model  of  the  structural  system.  In  general, 
accuracy  considerations  require  a  smaller  time  step  for  implicit  algorithms, 
as  compared  to  the  stability  requirements.  For  explicit  methods,  stability 
requirements  usually  dictate  the  maximum  time  step  size. 

41.  Bathe  and  Wilson  (1976)  describe  one  approach  used  to  assess  the 
accuracy  of  numerical  algorithm.  Illustrated  in  Figure  14  is  a  free  vibration 
problem  of  an  undamped  SDOF  system  subjected  to  an  initial  displacement  x0  and 
an  initial  velocity,  x  Figure  14a  shows  the  exact  solution  to  be  harmonic  with 
constant  maximum  amplitude.  Figure  14b  illustrates  results  for  the  free 
vibration  response  computed  using  a  numerical  algorithm.  The  errors  in  Figure 
14b  computed  results  take  the  form  of  amplitude  decay  with  time  and  period 
elongation  with  time,  as  described  in  Figure  14c.  The  order  of  accuracy  of 
numerical  algorithms  listed  in  Table  1  varies  in  proportion  to  the  square  of 
the  time  step  used  in  the  calculations. 

42.  One  approach  for  selecting  the  value  for  the  time  step  to  be  used 
in  the  calculations  is  to  require  a  prescribed  level  of  accuracy  for  the  com¬ 
puted  results.  This  decision  process  may  be  assisted  by  referring  to  charts 
developed  for  simplified  SDOF  system  problems  which  define  the  magnitude  of 
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(c)  Two  errors  associated  with  inaccurate  numerical 

solution 


Figure  14.  Example  of  inaccurate  response  for  an  undamped 
SDOF  system  in  free  vibration 
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the  errors  of  an  algorithm  (i.e.  amplitude  decay  and  period  elongation)  as  a 
function  of  the  time  step  size  At  used  in  the  calculation  and  the  natural 
period  T  of  the  numerical  model.  Examples  of  these  types  of  error  plots  for 
select  algorithms  are  shown  in  Figure  9.3  of  Bathe  and  Wilson  (1976)  and  in 
Figures  9.3.2  and  9.3.3  of  Hughes  (1987).  This  decision  for  a  complex  MDOF 
semidiscrete  finite  element  model  will  depend  upon  the  highest  (mode)  fre¬ 
quency  of  engineering  interest  for  the  structural  system  (see  Chapter  9  of 
Hughes  1987). 
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PART  VI:  RESPONSE  SPEC'"  JM  FOR  A  DAMPED  SDOF  SYSTEM 


43.  This  part  describes  the  construction  of  response  spectra  which  are 
graphs  of  the  maximum  values  of  acceleration,  velocity,  and/or  displacement 
response  of  an  infinite  series  of  damped  elastic  SDOF  systems  (Figure  10) 
subjected  to  an  acceleration  time  history  *ground(t).  These  maximum  response 
values  for  several  levels  of  damping  are  plotted  against  undamped  natural 
period  (units  of  seconds)  or  plotted  against  undamped  natural  cyclic  frequency 
of  vibration  (units  of  hertz  or  cycles/sec). 


Ground  Acceleration  Time  History 


44.  Figure  15  earthquake  accelerogram  will  be  used  as  the  ground  motion 
in  this  example.  It  is  the  South  00  degrees  East  (S00E)  horizontal  component 
recorded  at  the  El  Centro  site  in  Southern  California  during  the  Imperial 
Valley  Earthquake  of  18  May  1940  (Richter  magnitude  -  6.7).  The  recording 
station  is  founded  on  alluvium.  The  peak  ground  acceleration  is  equal  to 
341.7  cm/sec/sec  (0.35  g)  at  2.12  sec  into  ground  shaking.  The  ground 
velocity  and  displacement  curves  were  obtained  by  double  integration  of  the 
acceleration  time  history.  The  peak  velocity  and  peak  displacement  of  the 
ground  occur  at  2.19  sec  and  8.7  sec,  respectively.  It  is  typical  for  the 
three  peak  values  to  occur  at  different  times  during  earthquake  shaking. 

45.  Figure  15  accelerogram  started  out  as  an  analog  trace  recorded  by  a 
strong  motion  accelerograph .  The  acceleration  trace  vas  digitized  and  fil¬ 
tered  in  order  to  control  errors,  with  baseline  and  transducer  corrections 
applied  to  the  accelerogram.  Further  details  regarding  the  development  of  a 
corrected  accelerogram  are  described  in  Hudson  (1979).  The  California  Insti¬ 
tute  of  Technology  (Cal  Tech)  has  a  strong  motion  data  program  that  develops 
corrected  accelerograms  defined  at  0.02  sec  time  intervals.  The  corrected 
accelerograms  are  often  referred  to  as  the  Cal  Tech  Volume  II  corrected 
accelerograms . 


Response  Spectrum 


46.  A  response  spectrum  is  a  graphical  relationship  of  maximum  values 
of  acceleration,  velocity,  and/or  displacement  response  of  an  infinite  series 
of  linear  SDOF  systems  with  constant  damping  ratio  /?  shaken  by  the  same 
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IMPERIAL  VALLEY  EARTHQUAKE  MAY  18,  1940  -  2037  PST 

IIA001  40.001  .0  EL  CENTRO  SITE  IMPERIAL  VALLEY  IRRIGATION  DISTRICT  COMP  SOOE 
O  PEAK  VALUES  =  ACCEL  =  341.7  CM/SEC/SEt  VELOCITY  =  33.4  CM/SEC  DISPL  =  10.9  CM 
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Figure  15.  Ground  acceleration  and  integrated  ground  velocity  and 
displacement  time  histories,  1940  El  Centro  SOOE  component,  from 

Hudson  (1979) 


Figure  15  acceleration  time  history  in  this  example.  Each  Figure  10  SD0F  sys¬ 
tem  is  distinguished  by  the  value  selected  for  its  undamped  natural  period  of 
vibration  T  (Equation  13)  or  equivalently,  its  undamped  natural  cyclic  fre¬ 
quency  of  vibration  f  (Equation  14)  .  In  the  following  example  f)  is 
arbitrarily  set  equal  to  0.02. 


Peak  Response  Values  for  Each  SDOF  System 

47.  The  construction  of  the  response  spectrum  plots  a  succession  of 
SDOF  systems  with  fundamental  periods  T  ranging  from  near  zero  to  values  of 
several  seconds.  For  each  SDOF  system  of  value  T,  the  dynamic  response  is 
computed  using  one  of  the  numerical  procedures  listed  in  Table  1.  The  dynamic 
response  of  Figure  10  SDOF  system  is  expressed  in  terms  of  either  the  relative 
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response  or  the  total  response  of  the  SDOF  system.  Response  Spectrum  values 
are  the  maximum  response  values  for  each  of  the  five  types  of  SDOF  responses 
for  a  system  of  period  T  and  damping  0,  as  described  in  Table  2.  The  value 
assigned  to  each  of  the  five  dynamic  response  terms  for  a  SDOF  system  is  the 
peak  response  value  computed  during  earthquake  shaking. 


Relative  Displacement  Response  Spectrum 


48.  The  relative  displacement  response  spectrum,  SD  or  SD,  for  Fig¬ 
ure  15  ground  motion  is  shown  in  Figure  16,  computed  for  T  ranging  in  value 
from  0  sec  to  3  sec  and  using  0  -  0.02.  This  figure  displays  the  maximum 
absolute  relative  displacement  value,  |  x(t)  |max,  for  each  of  the  SDOF  sys¬ 
tems  analyzed,  as  illustrated  by  the  three  plots  showing  the  relative  dis¬ 
placement  time  histories  x(t)  that  were  computed  using  a  numerical  procedure 
for  the  SDOF  systems  having  T  equal  to  0.5  sec,  1.0  sec,  and  2.0  sec, 
respectively. 

SD  =  I  x(C)  |MX  (31) 


The  SD  values  for  the  three  SDOF  systems  are  listed  in  Table  3. 

49.  A  characteristic  of  relative  displacement  response  spectra  is  the 
absence  of  relative  movement  between  the  SDOF  system  and  the  ground  during 
earthquake  shaking  for  approximately  zero  values  of  T  (Figure  16) .  This  is 
due  to  the  fact  that  the  spring  stiffness  of  a  short  period,  high  frequency, 
SDOF  system  (refer  to  the  equation  for  T  in  Table  2)  is  so  large  that  it 
neither  stretches  nor  compresses  during  ground  shaking.  At  the  other  extreme, 
as  the  value  of  the  spring  stiffness  is  softened  (k  -*  0)  and  the  value  for  T 
is  greater  than  20  sec,  the  value  for  SD  approaches  the  maximum  ground  dis¬ 
placement  (not  shown  in  Figure  16) . 


Spectral  Pseudo-Velocifs 


50.  The  spectral  pseudo -veloc ity ,  Sv  or  PSV,  of  the  ground  motion 
fcgroundCt)  is  computed  using 


=  <*>  SD  =  SD 


38 


Table  2 

Definition  of  Earthquake  Response  Spectrum  Terms 


Symbols 

Definition 

DescriDtion 

SD  -  SD 

1  x(t)  1 

Relative  displacement  response  spectrum 
or  spectral  displacement 

SV 

1  x(t)  1  majt 

Relative  velocity  response  spectrum 

SA 

^  ^  max 

Absolute  acceleration  response  spectrum 

Sv  -  PSV 

c  2* 

Spectral  pseudo-velocity 

SA  -  PSA 

w2  s  -  ^  S 
w  bD  X2 

Spectral  pseudo-acceleration 

w  is  the  circular  frequency  of  vibration  of  the  undamped  SDOF  sysi ’u  in  units 
of  radians  per  second. 


T  is  the  natural  (or  fundamental)  period  of  vibration  of  the  undamped  SDOF 
system  in  units  of  seconds. 


for  each  of  the  SDOF  systems  analyzed.  Sv  is  distinguished  from  the  relative 
velocity  response  spectrum  SV,  to  be  described  in  a  subsequent  paragraph.  Sv 
for  Figure  15  ground  motion  shown  in  Figure  17b  was  computed  using  Equation  32 
for  each  T  ranging  in  value  from  0  to  3  sec  (at  closely  spaced  intervals)  in 
the  figure  (/9  -  0.02).  For  example,  with  SD  -  2.48  inches  for  the  T  -  0 . 5  sec 
(w  -  12.566  radians/sec)  SDOF  system,  Sv  -  31.16  inch/sec  by  Equation  32,  and 
identified  in  Figure  17b.  The  values  of  Sv  for  the  T  -  1.0  sec  and  T  - 
2.0  sec  SDOF  systems  are  listed  in  Table  3. 

51.  The  term  Sv  is  related  to  the  maximum  strain  energy  stored  within 
the  linear  spring  portion  of  the  SDOF  system  when  the  damping  force  is 
neglected , 
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EL  CENTRO, 
SOOE  COMPONENT, 


Figure  16.  Computation  of  displacement  response  spectrum 
f)  -  0.02,  1940  El  Centro  SOOE  component,  from 

Chopra  (1981) 
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Table  3 

Response  Spectral  Values  for  Three  SDOF  Systems 
with  8  —  0.02.  1940  El  Centro  S00E  Component 


T 

to 

SD 

Sy 

SA 

( sec) 

( radians/sec) 

( inch) 

( inch/sec) 

( inch/sec/sec) 

(R) 

0.5 

12.566 

2.48 

31.16 

391.62 

1.014 

1.0 

6.283 

6.61 

41.53 

260.95 

0.676 

2.0 

3.142 

8.84 

27.77 

87.25 

0.226 

w  is  the  circular  frequency  of  vibration  of  the  undamped  SDOF  system 


Ci)  = 


2n 

T 


where  T  is  the  natural  (or  fundamental)  period  of  vibration  of  the  undamped 
SDOF  system. 

The  spectral  pseudo-velocity .  Sv  -  w  SD 
The  spectral  pseudo -acceleration.  SA  -  w2  SD 
1  g  -  386.08858  inch/sec/sec 


Energy  =  -|  k  SD2  =  m  Sv2 


(33) 


Relative  Velocity  Response  Spectrum 

52.  The  relative  velocity  response  spectrum  SV  is  the  maximum  absolute 
value  of  the  computed  relative  velocity  time  history  (from  step  1  of  Fig¬ 
ure  10)  for  the  SDOF  system. 

SV  =  |  x  ( t)  |  (34) 

As  the  value  of  the  period  T  approaches  infinity  (i.e.  long  period  SDOF 
systems) ,  the  value  for  SV  approaches  the  maximum  ground  velocity 

I  Aground  (^)  I  max' 
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NATURAL  VIBRATION  PERIOD,  T,  sec 


(a) 


(b) 


(C) 


Figure  17.  Displacement,  pseudo-velocity 
pseudo-acceleration  linear  response 
spectra  plots,  /9  -  0.02,  1940  El  Centro 
S00E  component,  from  Chopra  (1981) 


Comparison  of  and  SV  Values 


53.  The  pseudo-velocity  value  Sv  for  a  SDOF  system  of  period  T  is  not 
equivalent  to  the  relative  velocity  value  SV,  as  shown  using  expressions 
derived  for  Sv  and  SV  in  Appendix  A  and  expressed  in  terms  of  the  Duharael's 
Integral.  However,  for  a  limited  range  of  periods,  the  value  for  Sv  is  an 
approximation  to  the  SV  value,  with  the  differences  between  the  values 
increasing  with  increasing  values  for  period  T.  Specifically,  for  short 
periods  (i.e.  high  frequencies)  the  values  for  Sv  are  nearly  the  same  as  the 
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values  for  SV.  For  long  periods  or  low  frequencies,  the  differences  in  the 
values  for  Sv  and  SV  can  be  significant. 

Spectral  Pseudo-Acceleration 


54.  The  spectral  pseudo-acceleration,  SA  or  PSA,  of  the  ground  motion 
xground ( t )  is  computed  using  either 


SA  =  <■>  Sy 


(35) 


or 


^  = 


or 


S*  = 


(36) 


for  each  of  the  SDOF  systems  analyzed.  SA  is  distinguished  from  the  absolute 
acceleration  response  spectrum  SA  to  be  described  in  a  subsequent  paragraph. 

SA  for  Figure  15  ground  motion,  shown  in  Figure  17c  was  computed  using  either 

Equation  35  or  Equation  36  for  each  T  ranging  in  value  from  0  to  3  sec  in  the 

figure  (/ 3  -  0.02%).  For  example,  with  SD  -  2.48  inches  for  the  T  -  0.5  sec 
(w  -  12.566  radians/sec)  SDOF  system,  SA  -  391.62  inch/sec/sec  (1.014  g)  by 
Equation  36  and  identified  in  Figure  17c.  The  values  of  SA  for  the  T  - 
1.0  sec  and  T  -  2.0  sec  SDOF  systems  are  listed  in  Table  3.  The  force  com¬ 
puted  as  the  product  of  the  mass  m  times  SA  is  a  good  approximation  to  the 

maximum  force  in  the  spring. 

Absolute  Acceleration  Response  Spectrum 

55.  The  absolute  acceleration  response  spectrum  SA  is  the  maximum  abso¬ 
lute  value  of  the  sum  of  the  computed  relative  acceleration  time  history  for 
the  SDOF  system  plus  the  ground  acceleration  time  history  (from  step  2  in 
Figure  10) . 

5A  =  |  xtotal(t)  |mx  =  j  x^^ft)  +*U)|mx  (37) 
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Comparison  of  S  and  SA  Values 

A  _ 


56.  The  pseudo -acceleration  value  SA  for  a  SDOF  system  of  period  T  is 
equal  to  the  absolute  acceleration  value  SA  when  /}  -  0  (refer  to  Appendix  A) . 
For  the  low  levels  of  damping  common  to  structural  dynamics  problems,  the 
values  for  SA  and  SA  are  nearly  equal  over  the  period  range  of  interest. 

57.  For  short  period,  high  frequency,  stiff  SDOF  systems  (T  <  0.05 
sec) ,  the  values  for  SA  and  SA  are  equal  to  the  peak  ground  acceleration 

I  Aground ( t )  |„,ax,  as  shown  in  Figure  18.  As  the  period  T  approaches  infinity 
(i.e.  massive  SDOF  systems),  the  values  for  SA  and  SA  approach  zero  (not  shown 
in  Figure  18) . 


Tripartite  Response  Spectra  Plot 

58.  A  tripartite  response  spectrum  plot  is  a  logarithmic  plot  of  the 
three  SDOF  response  quantities  SD,  Sv,  and  SA  as  a  function  of  '•he  undamped 
fundamental  period  T  in  units  of  seconds  (or  cyclic  frequency  f  in  units  of 
hertz) .  The  tripartite  response  spectrum  plot  for  Figure  15  ground  motion  is 
shown  in  Figure  19,  with  the  response  spectrum  values  computed  for  each  T 
ranging  in  value  from  0.04  to  3  sec  in  the  figure  (0  -  0.02).  The  plot  uses 
four  logarithmic  scales:  the  period  T  along  the  abscissa,  Sv  along  the 
ordinate,  SA  along  an  axis  oriented  at  45  deg  counterclockwise  from  hori¬ 
zontal,  and  SD  along  an  axis  oriented  at  45  deg  clockwise  from  horizontal. 

The  response  quantities  SD  -  2.48  inches,  Sv  -  31.16  inch/sec,  and  SA  - 
1.014  g  for  the  T  -  0.5  sec  SDOF  system  identified  in  Figure  19. 

59.  Figure  19  tripartite  response  spectrum  plot  has  the  advantage  of 
condensing  the  information  presented  using  three  plots  in  Figure  17  onto  a 
single  plot  because  of  the  interrelationships  between  the  three  terms  (refer 
to  Equations  32,  35,  and  36). 

60.  The  tripartite  response  spectrum  plot  for  the  1940  El  Centro  ground 
motion  (Figure  15)  in  Figure  20  shows  the  response  spectrum  values  for  levels 
of  damping,  0-0,  0.02,  0.05,  0.1  and  0.2.  This  figure  illustrates  the  trend 
common  to  all  ground  motions  of  increased  values  for  SD,  Sv,  and  SA  with  lower 
damping  levels. 


Figure  18.  Pseudo-acceleration  and  absolute  acceleration  linear  response 

spectrum  plots 


Fundamental  Periods  at  Which  the  Maximum  Values  for 


S  ,  S  ,  and  S  are  Computed 

A  V  U 


61.  For  the  1940  El  Centro  ground  motion,  the  largest  value  of  SA  equal 
to  1.29  g  (1,265  cm/sec/sec),  for  an  SDOF  system  with  T  equal  to  0.47  sec  (f  - 
2.13  Hz)  and  0  -  0.02  (refer  to  Figure  17  or  Figures  19  and  20).  This  value 
is  3.7  times  larger  than  the  peak  acceleration  value  of  0.35  g  (341.7  cm/sec/ 
sec)  of  the  accelerogram  (Figure  15).  The  largest  value  for  Sv,  equal  to 
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Figure  19.  Tripartite  logarithmic  plot  of  response  spectrum, 
$  -  0.02,  1940  EL  Centro  SOOE  component,  from  Chopra  (1981) 
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RESPONSE  SPECTRUM 

IMPERIAL  VALLEY  EARTHQUAKE  MAY  18.  1940  -  2037  PST 
IIIA001  40.001.0  EL  CENTRO  SITE  IMPERIAL  VALLEY  IRRIGATION  DISTRICT  COMP  SOOE 
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Figure  20.  Tripartite  logarithmic  plot  of  response  spectrum,  0  **  0 , 
0.02,  0.05,  0.1,  and  0.2,  1940  El  Centro  SOOE  Component,  from 

Chopra  (1981) 


4 5.7  in. /sec  (18  cm/sec),  is  computed  for  a  value  of  T  equal  to  0.84  sec  (f  - 
1.19  Hz).  The  largest  value  for  SD,  equal  to  18  in.  (7.09  cm),  is  computed 
for  a  value  of  T  equal  to  10.5  sec  (f  -  0.095  Hz)  (Figure  20).  The  order  of 
increasing  T  values  (seconds)  with  the  peak  values  of  SA,  Sv,  and  SD  is  typi¬ 
cal  of  response  spectra  for  ground  motions.  The  order  of  spectral  values  is 
reversed  when  the  SDOF  system  is  described  in  terms  of  increasing  values  of 
cyclic  frequency  f  (Hz)  (Equation  14). 
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Frequency  Regions  of  the  Tripartite  Response  Spectra 


62.  Although  quite  irregular  with  local  peaks  and  spikes,  the  1940  El 
Centro  earthquake  motion  tripartite  response  spectrum  has  the  general  shape  of 
a  trapezoid  or  tent  (refer  to  Figure  20).  This  general  shape  implies  three 
distinct  regions  with  regards  to  the  frequency  content  of  the  1940  El  Centro 
earthquake  motion. 

63.  Within  the  short  period  range  of  1/8  sec  <  T  <  1/2  sec  or 
equivalently,  the  high  frequency  range  of  2  Hz  <  f  <  8  Hz,  SA  is  nearly  con¬ 
stant  (an  average  value  of  0.95  g  (931.6  cm/sec/sec)  with  a  variation  equal  to 
±  25%)  ,  as  compared  to  the  variation  in  Sv  and  SD  values  (Figure  19)  . 

64.  There  are  two  distinct  subregions  within  the  intermediate  period 
range  of  the  1940  El  Centro  motion  response  spectra  shown  in  Figure  19.  The 
first  subregion  is  within  the  period  range  of  1/2  sec  <  T  <  1  sec  (1  Hz  <  f  < 

2  Hz) ,  with  Sv  a  nearly  constant  38  in/sec  (15  cm/sec)  value  (±  20%) ,  as  com¬ 
pared  to  the  variation  in  SA  and  SD  values.  Between  1  sec  <  T  <  3  sec  (or 
0.33  Hz  <  f  <  1  Hz),  Sv  is  nearly  a  constant  28  in. /sec  (11  cm/sec)  value 

(±  25%). 

65.  The  long  period  range,  T  greater  than  3  sec  (the  low  frequency 
range  of  f  <  0.33  Hz),  SD  averages  13  in  (5.1  cm)  for  0  -  0.02,  with  a  varia¬ 
tion  equal  to  ±  40%  (Figure  20).  For  this  accelerogram,  the  variation  in  the 
average  value  of  SD  is  greater  than  the  variations  in  the  SA  and  Sv  values  for 
the  short  and  intermediate  period  ranges,  respectively. 

Design  Response  Spectra 

66.  The  previously  described  trends  in  the  spectral  content  with  fre¬ 
quency  of  earthquake  accelerograms  were  recognized  in  the  1970 's  by  numerous 
earthquake  engineers  and  seismologists.  Statistical  analyses  of  response 
spectra  of  earthquake  motions  were  conducted  by  several  groups  with  the 
objectives  of  identifying  the  factors  affecting  the  response  spectra  for 
earthquake  motions,  characterizing  shape  or  the  frequency  content  of  the 
earthquake  spectrum  for  the  category  under  study  and  developing  smooth,  broad 
band  spectrum  for  use  in  the  design  of  structures  for  earthquake  loadings.  A 
broad  band  spectrum  ensures  that  sufficient  seismic  energy  is  delivered  to  all 
frequencies . 
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67.  The  results  of  some  of  the  early  ground  motion  studies  are  summa¬ 
rized  in  Seed,  Ugas ,  and  Lysmer  (1976),  Mohraz  (1976),  and  Newmark  and  Hall 
(1982).  Their  results  show  that  the  spectra  frequency  content  of  the  recorded 
accelerograms  are  dependent  upon  the  earthquake  magnitude,  distance  from 
causative  fault  to  site,  the  type  of  substratum  on  which  the  recording  station 
is  founded  (i.e.,  rock,  shallow  alluvium,  or  deep  alluvium),  and  the  tectonic 
environment.  Using  standard  regression  techniques  for  selected  ground  motion 
categories,  (smooth)  trends  in  the  variation  in  spectral  content  with 
frequency  were  quantitatively  identified,  as  well  as  their  variation  about 
their  median  values,  expressed  in  terms  of  a  standard  deviation. 


Newmark  and  Hall  Design  Response  Spectrum 


68.  The  Newmark  and  Hall  (1982)  design  response  spectrum  exemplifies 
the  important  features  common  to  current  procedures  used  to  characterize  the 
spectral  content  of  earthquake  motions.  A  Newmark  and  Hall  tripartite  design 
response  spectrum  for  a  moderate  earthquake  at  a  competent  soil  site  corre¬ 
sponds  to  the  solid  line  shown  in  Figure  21  for  -  0.05.  The  peak  ground 
acceleration,  peak  ground  velocity,  and  peak  ground  displacement  used  in  this 
example  and  included  for  reference  in  this  figure  (dashed  line),  are  0.5  g 
(490.3  cm/sec/sec),  61  cm/sec,  and  45  cm,  respectively.  The  spectral  values 
are  plotted  as  a  function  of  the  frequency  f  (Hz),  rather  than  period  T  as 
Figures  19  and  20.  The  Newmark  and  Hall  design  response  spectrum  is  a  smooth, 
broad  band  spectrum  composed  of  three  distinct  regions  of  constant,  limiting 
spectral  values  for  SA,  Sv,  and  SD,  as  identified  in  the  figure.  These  spec¬ 
tral  amplitude  values  correspond  to  one  standard  deviation  above  the  median 
spectral  values  (84.1  %  cumulative  probability).  Mean  spectral  values  are 
presented  in  the  Newmark  and  Hall  (1982)  publication  as  well.  A  smooth,  broad 
band  spectrum  ensures  that  sufficient  seismic  energy  is  delivered  to  all  fre¬ 
quencies  in  a  dynamic  analysis  of  a  MDOF  semidiscrete  structural  model.  For 
details  regarding  the  construction  of  the  Newmark  and  Hall  design  response 
spectrum,  consult  their  1982  EERI  monograph. 


Response  Spectrum  in  the  Dynamic  Analysis  of  MDOF  Structural  Models 


69.  The  response  spectrum  is  used  to  represent  the  seismic  loading  in 
the  dynamic  analysis  of  MDOF  semidiscrete  structural  models.  For  details 
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FREQUENCY,  Hz 

Figure  21.  Newmark  and  Hall  elastic  design  response 
spectrum,  from  Newmark  and  Hall  (1982) 

regarding  its  application  in  structural  dynamic  problems,  the  reader  is 
referred  to  textbooks  such  as  Clough  and  Penzien  (1975),  Chopra  (1981),  or  Paz 
(1985)  . 
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APPENDIX  A:  SOLUTIONS  FOR  SPECTRAL  TERMS  BASED  UPON  DUHAMEL'S 
INTEGRAL  SOLUTION  FOR  RELATIVE  DISPLACEMENT 

1.  This  appendix  describes  the  derivation  of  relationships  for  each  of 
the  five  response  spectrum  terms  SD,  Sv,  SV,  SA>  and  SA  (Table  2),  expressed 
in  terms  of  Duhamel's  integral.  The  resulting  relationships  show  that  the 
terms  Sv  and  SV  are  not  equivalent,  while  the  terms  SA  and  SA  are  equivalent 
only  if  0  -  0. 

Duhamel's  Integral  Solution  for  Relative  Displacement 

2.  The  relative  displacement  of  the  SDOF  system  was  described  in  Part 
IV  involving  Duhamel's  integral.  The  derivation  of  x(t)  (Equation  22) 
involved  the  representation  of  the  load  time  history  P(t)  -  -  mxRround  as  a 
series  of  loadings  P(r)-  -  mxground  applied  to  the  SDOF  system  for  short  dura¬ 
tion  time  intervals  dr .  The  resulting  relationship  for  a  damped  SDOF  system 
is  given  as 

x(t)  =  -  -L  f'  xground(x)  e-*w‘e_t>  sin  [wD(t  -  x)  ]  dx  (Al) 


where 

wd  -  the  damped  angular  frequency  of  vibration  (Equation  16) 
w  -  the  undamped  angular  frequency  of  vibration  (Equation  11) 
P  =  the  fraction  of  critical  damping  (Equation  12) . 


Relative  Displacement  Response  Spectrum  S^ 
Expressed  in  Terms  of  Duhamel's  Integral 


3.  The  relative  displacement  response  spectrum  SD  is  equal  to  the  maxi¬ 
mum  absolute  relative  displacement  value  |  x(t)  |max  (Equation  25)  for  a 
damped  SDOF  system  of  period  T.  Introducing  Equation  Al  for  the  relati\ a 
displacement  x(t),  SD  becomes 


= 


-  -  r 


*“  ground 


(x)  sin  [o>D(t  -  x)]  dx 


(A2) 


Al 


Spectral  Pseudo-Velocity  Expressed  In  Terms  of  Duhamel's  Integral 


4.  Sv  may  be  expressed  in  terms  of  Duhamel's  integral  by  introducing 
the  Duhamel's  integral  solution  for  SD  into  Equation  26,  resulting  in  the 
relationship  for  Sv  is  given  by 

Sv  =  ~  7T  /'  Sgioimd^)  e->“(e-t)  sin  [*»B(t-t)]  eft  (A! 


For  an  undamped  SDOF  system  (0-0),  this  expression  simplifies  to 


=  I  -  jC  *ground<T>  Sin  t«  (  t  -  t)]  eft 


Relative  Velocity  Response  Spectrum  SV  Expressed 
in  Terms  of  Duhamel's  Integral 

5.  The  relative  velocity  response  spectrum  SV  is  equal  to  the  maximum 
absolute  relative  velocity  value  |  x(t)  |max  (Equation  28)  for  a  damped  SDOF 
system  of  period  T.  Differentiating  Duhamel's  integral  solution  for  relative 
displacement  (Equation  Al)  with  respect  to  time  results  in  an  expression  for 
x(t),  with  SV  equal  to  the  maximum  absolute  value  of  the  resulting  expression. 


ground**)  e'""'*  '  *’  COS  [uc(t  -t)]  dx 


,(t)  sin  [ut(t  -  t)]  efi 


Refer  to  page  62  of  Hudson  (1979)  or  section  1.5  of  Newmark  and  Rosenbluth 
(1971)*.  For  an  undamped  SDOF  system  (0  -  0)  Equation  A5  becomes 


*3round<T>  COS  [«(t  -  T)]  dx 


*  References  cited  in  this  appendix  are  included  in  the  References  at  the  end 
of  the  main  text. 
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Comparison  of  Spectral  Pseudo-Velocitv  and  Relative  Velocity 
Response  Spectrum  Relationships 

6.  Comparison  of  the  expressions  for  Sv  (Equation  A3)  and  SV 
(Equation  A5)  in  the  previous  two  paragraphs  demonstrates  the  difference 
between  the  two  terms.  The  simplified  relationships  for  an  undamped  SDOF 
system  (/9  -  0)  shows  the  terms  for  Sv  (Equation  A4)  and  SV  (Equation  A6)  dif¬ 
fer  by  the  trigonometric  function  used  within  Duhamel’s  integral. 

Spectral  Pseudo-Acceleration  SA  Expressed  in 
Terms  of  Duhamel's  Integral 


7.  SA  may  be  expressed  in  terms  of  Duhamel's  integral  by  introducing 
the  Duhamel's  integral  solution  for  SD  into  Equation  30,  resulting  in  the 
relationship  for  SA  given  by 


SA  = 


-  /C*ground(t)  Sin  [<*>p ( c  —  t )  3  ctx 


(A7) 


For  an  undamped  SDOF  system  (0  -  0) ,  this  expression  simplifies  to 


SA  = 


-  w  /"  xgtound(T)  sin  [«(t  -  t)]  dx 


(A8) 


Absolute  Acceleration  Response  Spectrum  in 
Terms  of  Duhamel's  Integral 

8.  The  absolute  acceleration  response  spectrum  SA  is  the  maximum 
absolute  value  of  the  sum  of  the  computed  relative  acceleration  time  history 
for  the  SDOF  system  plus  the  ground  acceleration  time  history. 

SA  =  |  xtotal(t)  |mx  =  j  xgroundU)  +x(t)|max  (A9) 

To  introduce  Duhamel's  integral  solution  for  relative  displacement  into  this 
relationship,  Xtotai^)  is  expressed  in  terms  of  the  relative  velocity  and 
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relative  displacement  of  the  SDOF  system.  Rearranging  the  dynamic  equilibrium 
Equation  21  and  introducing 

*to«l  =  abound  *  *  <A10> 


results  in  the  relationship 

=  -2P«x-<o2JC 


(All) 


Introducing  Equation  Al  and  the  term  within  the  absolute  value  in  Equation  A5 
into  Equation  All,  results  in  the  relationship 


2wp  ^ground  e'p“<£  *  cos  [«fl(t-x)]  dx 
0 

/C*gr0UIld( X)  sin  [»D(t-t)]  dil 

/I-T  o  1 


(A12) 


Refer  to  page  63  of  Hudson  (1979)  or  section  1.5  of  Newmark  and  Rosenbluth 
(1971).  For  an  undamped  SDOF  system  (f)  -  0)  Equation  A12  becomes 


SA  = 


•/* 


Aground 


( x)  sin  [» ( t  -  t)  )  dx 


(A13) 


Comparison  of  Spectral  Pseudo-Acceleration  and  Absolute 
Acceleration  Response  Spectrum  Relationships 

9.  Comparison  of  the  expressions  for  SA  (Equation  A7)  and  SA  (Equa¬ 
tion  A12)  in  the  previous  two  paragraphs  demonstrates  the  difference  between 
the  two  terms.  However,  the  simplified  relationships  for  an  undamped  SDOF 
system  ( fi  —  0)  for  SA  (Equation  A8)  and  SA  (Equation  A13)  are  the  same. 
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Theoretical  Basis  for  CTABS80:  A  Computer  Program  for 
Three-Dimensional  Analysis  of  Building  Systems 

User's  Guide:  Computer  Program  for  Analysis  of  Beam-Column 
Structures  with  Nonlinear  Supports  (CBEAMC) 


Feb  1978 
Mar  1979 

Jan  1980 
Jan  1 980 

Feb  1980 

Mar  1980 

Jun  1980 
Jun  1982 
Aug  1983 

Dec  1980 
Dec  1980 

Dec  1 980 
Dec  1 980 

Dec  1 980 


Feb  1981 
Mar  1981 

Feb  1981 
Mar  1981 
Mar  1981 

Mar  1981 
Aug  1981 
Sep  1981 
Jun  1982 
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Instruction  Report  K-82-7 
Instruction  Report  K-83-1 
Instruction  Report  K-83-2 
Instruction  Report  K-83-5 

Technical  Report  K-83-1 
Technical  Report  K-83-3 

Technical  Report  K-83-4 
Instruction  Report  K-84-2 

Instruction  Report  K-84-7 

Instruction  Report  K-84-8 

Instruction  Report  K-84-1 1 

Technical  Report  K-84-3 

Technical  Report  ATC-86-5 

Technical  Report  ITL-87-2 
Instruction  Report  ITL-87-1 
Instruction  Report  ITL-87-2 
Technical  Report  ITL-87-6 
Instruction  Report  ITL-87-3 


_ —  (Continued) 

Title 

User’s  Guide:  Computer  Program  for  Bearing  Capacity  Analysis 
of  Shallow  Foundations  (CBEAR) 

User’s  Guide:  Computer  Program  with  Interactive  Graphics  for 
Analysis  of  Plane  Frame  Structures  (CFRAME) 

User’s  Guide:  Computer  Program  for  Generation  of  Engineering 
Geometry  (SKETCH) 

User's  Guide:  Computer  Program  to  Calculate  Shear,  Moment, 
and  Thrust  (CSMT)  from  Stress  Results  of  a  Two-Dimensional 
Finite  Element  Analysis 

Basic  Pile  Group  Behavior 

Reference  Manual:  Computer  Graphics  Program  for  Generation  of 
Engineering  Geometry  (SKETCH) 

Case  Study  of  Six  Major  General-Purpose  Finite  Element  Programs 

User’s  Guide:  Computer  Program  for  Optimum  Dynamic  Design 
of  Nonlinear  Metal  Plates  Under  Blast  Loading  (CSDOOR) 

User's  Guide:  Computer  Program  for  Determining  Induced 
Stresses  and  Consolidation  Settlements  (CSETT) 

Seepage  Analysis  of  Confined  Flow  Problems  by  the  Method  of 
Fragments  (CFRAG) 

User’s  Guide  for  Computer  Program  CGFAG,  Concrete  General 
Flexure  Analysis  with  Graphics 

Computer-Aided  Drafting  and  Design  for  Corps  Structural 
Engineers 

Decision  Logic  Table  Formulation  of  ACI  318-77,  Building  Code 
Requirements  for  Reinforced  Concrete  for  Automated  Con¬ 
straint  Processing,  Volumes  I  and  II 

A  Case  Committee  Study  of  Finite  Element  Analysis  of  Concrete 
Flat  Slabs 

User’s  Guide:  Computer  Program  for  Two-Dimensional  Analysis 
of  U-Frame  Structures  (CUFRAM) 

User’s  Guide:  For  Concrete  Strength  Investigation  and  Design 
(CASTR)  in  Accordance  with  ACI  318-83 

Finite-Element  Method  Package  for  Solving  Steady-State  Seepage 
Problems 

User's  Guide:  A  Three  Dimensional  Stability  Analysis/Design 
Program  (3DSAD)  Module 
Report  1 :  Revision  1 :  General  Geometry 
Report  2:  General  Loads  Module 
Report  6:  Free-Body  Module 


Date 

Jun  1982 
Jan  1 983 
Jun  1983 

Jul  1983 

Sep  1983 
Sep  1 983 

Oct  1983 
Jan  1984 

Aug  1984 

Sep  1984 

Sep  1984 

Oct  1984 

Jun  1986 

Jan  1987 
Apr  1987 
May  1987 
May  1987 

Jun  1987 

Jun  1987 
Sep  1989 
Sep  1989 
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Title 


Date 


Instruction  Report  ITL-87-4 
Technical  Report  ITL-87-4 


Instruction  Report  GL-87-1 

Instruction  Report  ITL-87-5 
Instruction  Report  ITL-87-6 

Technical  Report  ITL-87-8 

Instruction  Report  ITL-88-1 

Technical  Report  ITL-88-1 

Technical  Report  ITL-88-2 

Instruction  Report  ITL-88-2 

Instruction  Report  ITL-88-4 

Instruction  Report  GL-87-1 

Technical  Report  ITL-89-3 
Technical  Report  ITL-89-4 


User's  Guide:  2-D  Frame  Analysis  Link  Program  (LINK2D) 


Finite  Element  Studies  of  a  Horizontally  Framed  Miter  Gate 

Report  1 :  Initial  and  Refined  Finite  Element  Models  (Phases 
A,  B,  and  C),  Volumes  I  and  II 
Report  2:  Simplified  Frame  Model  (Phase  D) 

Report  3:  Alternate  Configuration  Miter  Gate  Finite  Element 
Studies-Open  Section 

Report  4:  Alternate  Configuration  Miter  Gate  Finite  Element 
Studies-Closed  Sections 

Report  5:  Alternate  Configuration  Miter  Gate  Finite  Element 
Studies-Additiona!  Closed  Sections 
Report  6:  Elastic  Buckling  of  Girders  in  Horizontally  Framed 
Miter  Gates 

Report  7:  Application  and  Summary 


User's  Guide:  UTEXAS2  Slope-Stability  Package:  Volume  I, 
User’s  Manual 


Sliding  Stability  of  Concrete  Structures  (CSLIDE) 

Criteria  Specifications  for  and  Validation  of  a  Computer  Program 
for  the  Design  or  Investigation  of  Horizontally  Framed  Miter 
Gates  (CMITER) 

Procedure  for  Static  Analysis  of  Gravity  Dams  Using  the  Finite 
Element  Method  -  Phase  la 

User’s  Guide:  Computer  Program  for  Analysis  of  Planar  Grid 
Structures  (CGRID) 

Development  of  Design  Formulas  for  Ribbed  Mat  Foundations 
on  Expansive  Soils 

User’s  Guide:  Pile  Group  Graphics  Display  (CPGG)  Post¬ 
processor  to  CPGA  Program 

User's  Guide  for  Design  and  Investigation  of  Horizontally  Framed 
Miter  Gates  (CMITER) 


User's  Guide  for  Revised  Computer  Program  to  Calculate  Shear, 
Moment,  and  Thrust  (CSMT) 

User’s  Guide:  UTEXAS2  Slope-Stability  Package;  Volume  II, 
Theory 

User’s  Guide:  Pile  Group  Analysis  (CPGA)  Computer  Group 

CBASIN-Structural  Design  of  Saint  Anthony  Falls  Stilling  Basins 
According  to  Corps  of  Engineers  Criteria  for  Hydraulic 
Structures;  Computer  Program  X0098 


Jun  1987 
Aug  1987 


Aug  1987 

Oct  1987 
Dec  1987 

Jan  1988 

Feb  1988 

Apr  1988 

Apr  1988 

Jun  1988 

Sep  1988 

Feb  1989 

Jut  1989 
Aug  1989 
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Technical  Report  ITL-89-5 

Technical  Report  ITL-89-6 
Contract  Report  ITL  89-1 
Instruction  Report  ITL-90-1 
Technical  Report  ITL-90-3 

Instruction  Report  ITL-90-6 
Instruction  Report  ITl-90-2 
Technical  Report  ITL-91-3 

Instruction  Report  ITL-91-1 

Instruction  Report  ITL-87-2 
(Revised) 

Technical  Report  ITL-92-2 
Technical  Report  ITL-92-4 


(Concluded) 

Title 

CCHAN-Structural  Design  of  Rectangular  Channels  According 
to  Corps  of  Engineers  Criteria  for  Hydraulic 
Structures;  Computer  Program  X0097 

The  Response-Spectrum  Dynamic  Analysis  of  Gravity  Dams  Using 
the  Finite  Element  Method;  Phase  II 

State  of  the  Art  on  Expert  Systems  Applications  in  Design, 
Construction,  and  Maintenance  of  Structures 

User's  Guide;  Computer  Program  for  Design  and  Analysis 
of  Sheet  Pile  Walls  by  Classical  Methods  (CWALSHT) 

Investigation  and  Design  of  U-Frame  Structures  Using 
Program  CUFRBC 

Volume  A:  Program  Critena  and  Documentation 
Volume  B:  User's  Guide  for  Basins 
Volume  C:  User's  Guide  for  Channels 

User's  Guide:  Computer  Program  for  Two-Dimensional  Analysis 
of  U-Frame  or  W-Frame  Structures  (CWFRAM) 

User's  Guide:  Pile  Group-Concrete  Pile  Analysis  Program 
(CPGC)  Preprocessor  to  CPGA  Program 

Application  of  Finite  Element,  Grid  Generation,  and  Scientific 
Visualization  Techniques  to  2-D  and  3-D  Seepage  and 
Groundwater  Modeling 

User's  Guide.  Computer  Program  for  Design  and  Analysis 
of  Sheet-Pile  Walls  by  Classical  Methods  (CWALSHT) 

Including  Rowe's  Moment  Reduction 

User's  Guide  for  Concrete  Strength  Investigation  and  Design 
(CASTR)  m  Accordance  with  ACI  318-89 

Fiimte  Element  Modeling  of  Welded  Thick  Plates  for  Bonneville 
Navigation  Lock 

Introduction  to  the  Computation  of  Response  Spectrum  for 
Earthquake  Loading 


Date 

Aug  1989 

Aug  1989 
Sep  1989 
Feb  1990 
May  1990 

Sep  1 990 
Jun  1990 
Sep  1990 

Oct  1991 

Mar  1992 
May  1992 
June  1992 


